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The focus of this article is on recent progress in muon spin rotation (/iSR) studies of the vortex 
cores in type-II superconductors. By comparison of itSR measurements of the vortex-core size in a 
variety of materials with results from techniques that directly probe electronic states, the effect of 
delocalized quasiparticles on the spatial variation of field in a lattice of interacting vortices has been 
determined for both single-band and multi-band superconductors. These studies demonstrate the 
remarkable accuracy of what some still consider an exotic technique. In recent years /iSR has also 
been used to search for magnetism in and around the vortex cores of high-temperature superconduc- 
tors. As a local probe pSR is specially suited for detecting static or quasistatic magnetism having 
short-range or random spatial correlations. As discussed in this review, iiSR experiments support a 
generic phase diagram of competing superconducting and magnetic order parameters, characterized 
by a quantum phase transition to a state where the competing order is spatially non-uniform. 



INTRODUCTION 

Muon spin rotation (/iSR) is an experimental tech- 
nique primarily used to measure local magnetic fields 
inside samples. The discovery of high-transition temper- 
ature (high-T c ) superconductivity in 1986 brought about 
a rapid world-wide expansion in the use and applications 
of ^iSR. Since then zxSR has been routinely applied to to 
investigations of these and other newly discovered type-II 
superconductors. The technique allows for studies in zero 
external magnetic field, which combined with its sensitiv- 
ity as a local probe has provided distinctive information 
on the occurrence of internal magnetism as a co-existing 
or competing phase, or as a consequence of time-reversal 
symmetry breaking superconductivity. From zero-field 
/iSR studies of cuprates, a generic temperature- vs-doping 
phase diagram has been constructed, showing the coexis- 
tence of high-T c superconductivity with static magnetism 
in lightly doped samples. Today there is still much de- 
bate on the origin of this magnetism and its importance 
to the high-T c 'problem'. 

The vortex state provides another avenue for investiga- 
tion of type-II superconductors with /iSR [l| . For many 
years such studies focussed solely on obtaining experi- 
mental information on the magnetic penetration depth 
(A), through measurements of the muon spin depolariza- 
tion rate (a) resulting from the broad internal magnetic 
field distribution n{B) of the flux-line lattice (FLL). The 
temperature and magnetic field dependences of A, which 
in many systems can also be determined in the Meiss- 
ner phase by other techniques, reflect the pairing state 
symmetry of the superconducting carriers. With further 
advances of the /iSR method came the ability to focus at- 
tention on the properties of the vortex cores themselves. 

The first-ever study to account for the finite size of the 
vortex cores in the analysis of /xSR data was an inves- 
tigation by Herlach et al. of n(B) in pure Nb single 
crystals. The measured field distributions were shown to 



be consistent with numerical solutions of the microscopic 
BCS-Gor'kov theory. Some years later, the magnetic field 
dependence of the vortex core size was determined from 
ixSR measurements on single-crystal NbSe2 H. The re- 
sults confirmed earlier scanning tunnelling spectroscopy 
(STS) measurements on NbSe2 that showed a shrinking 
of the vortex cores with increasing magnetic field jij. 
This behaviour could be attributed to an increased over- 
lap of the quasiparticle states around a vortex core with 
those coming from neighbouring vortices. However, these 
/iSR studies were more than just another means of access- 
ing information obtainable by another experimental tech- 
nique. Instead they marked the development of a more 
powerful method for investigating some of the intrinsic 
properties of vortex cores in type-II superconductors. 

The STS technique, which is sensitive to the elec- 
tronic structure of the vortex cores, is limited to probing 
individual vortices near the sample surface. Near the 
surface the vortices spread out [H, H| and their proper- 
ties are strongly influenced by surface inhomogeneities 
and/or defects. Today one can study vortices immedi- 
ately above or below the surface by i/SR using low-energy 
(several keV) positively charged muons [1,0], or by 
/3-detected NMR using low-energy radioactive ions [8(]. 
In contrast, the experiments of Refs. 0, [j| used ener- 
getic (~ 4 MeV) t< + that stop at interstitial or bond sites 
in the bulk of the sample where they directly probe the 
local magnetic fields. The term "bulk" means that the 
stopping range of these faster muons is approximately 
150 mg/cm 2 , which requires samples ~ 1 mm thick. In 
further contrast to the STS method, /iSR studies yield 
average information on the vortex cores, using ~ 10 7 /z + 
to randomly probe the ~ 10 9 vortices in a typical size 
sample. 

Since the experiments of Ref. [|[ , strong field and tem- 
perature dependences of the vortex-core size have been 
found by /iSR in a variety of superconductors. Through 
comparison of the results with theoretical models and 
experiments that are directly sensitive to quasiparticles 
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properties, a good understanding of many of the fj.SH ex- 
periments has been achieved. In s-wave superconductors, 
it is now well established that the vortex core size de- 
pends on both the thermal occupancy of the bound quasi- 
particle core states and the overlap of the correspond- 
ing quasiparticle wave functions with those of nearest- 
neighbour vortices. However, in exotic systems such 
as high-T c superconductors, where localized cores states 
may be absent, there is currently insufficient experimen- 
tal information to make similar definitive statements. On 
the other hand, recent fj.SH studies of the vortex cores in 
underdoped high-T c superconductors have shed new light 
on the ground state that emerges when superconductiv- 
ity is suppressed. Combining information obtained from 
/iSR experiments in zero and nonzero magnetic fields, 
the latest results support a picture of closely competing 
superconducting and magnetic ground states. 
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MUON SPIN ROTATION (fj.SH) 

General description 

The acronym fiSK dates back to 1974, and stands for 
either 'muon spin relaxation, rotation' or 'resonance'. 
These three terms describe different uses of the magnetic 
moment of a muon or fi + ) to probe matter. While 
the principle aspects of the technique are analogous to 
nuclear magnetic resonance (NMR), there are many im- 
portant differences. Only some of these are touched upon 
in this review. 

The primary use and strength of /iSR is its unmatched 
sensitivity to internal magnetism. Central to the /zSR 
method is the use of a nearly 100 % spin-polarized muon 
beam, naturally generated from the weak interaction de- 
cay of pions. This is a great advantage over conventional 
NMR, which relies on thermal equilibrium nuclear spin 
polarization in a large magnetic field. Zero-field (ZF) 
/iSR is routinely used to study small internal magnetic 
fields of natural origin. In contrast to neutron scat- 
tering, the information provided by /xSR is integrated 
over reciprocal space, which makes it ideal for studies 
of short-range magnetic correlations or disordered mag- 
netism. The magnetic moment of the muon is 3.18 times 
larger than that of the proton, making it even more sensi- 
tive to magnetism than NMR. Although generally a nui- 
sance in experiments, /iSR even detects the dipolar fields 
of nuclear moments. In fact, magnetic fields as small 
as ~ 0.1 G are detectable — although it is important to 
emphasize that this refers to the local field at the muon 
stopping site. 



FIG. 1: Schematic of the arrangement for a TF-/^SR exper- 
iment on a type-II superconductor in the vortex state. The 
muon spin Larmor precesses about the local magnetic field 
B at its stopping site in the sample, and subsequently un- 
dergoes the three-body decay fj + — > e + + v e + The time 
evolution of the muon spin polarization P x (t) is accurately 
determined by detection of the decay positrons from ~ 10 6 
muons implanted one at a time. 

Transverse-field fiSH. 

The internal magnetic field distribution n(B) of a type- 
II superconductor in the vortex state is measured by the 
so-called 'transverse-field' muon spin rotation (TF-/iSR) 
method. The geometry of a TF-/iSR experiment is shown 
in figure [TJ The external magnetic field H is applied 
transverse to the direction of the initial muon spin polar- 
ization P x (0) 7 which defines the x-axis. In high-T c super- 
conductors, the positive muon (^ + ) forms an ~ 1 A bond 
with an oxygen atom [^, fioj ] . but in general the muon 
will stop at an interstitial site in the sample. There the 
muon spin precesses about the local magnetic field B(r) 
in a plane perpendicular to the local field axis. The muon 
subsequently decays, emitting a fast positron. The an- 
gular dependence of the decay probability of the muon is 
given by 

W(E, 6) = 1 + a{E) cos(0) , (1) 

where E is the kinetic energy of the decay positron, 6 is 
the angle between the directions of the muon spin and 
the emitted positron, and a(E) is an asymmetry factor. 
When all positron energies are sampled with equal prob- 
ability, the asymmetry factor has the value a = 1/3. The 
statistical average direction of the muon spin polarization 
is obtained by measuring the anisotropic angular distri- 
bution of decay positrons from an ensemble of implanted 
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nations. 

The /iSR signal obtained by the detection of the decay 
positrons is given by 

A(t) = a P x (t) , (2) 

where A(t) is the /iSR 'asymmetry' spectrum, ao is the 
asymmetry maximum, and P x (t) is the time evolution of 
the muon spin polarization 

P x (t) = / n(B) cos^Bt + <P)dB . (3) 
Jo 

Here 7 P = 0.0852 /is -1 G _1 is the muon gyromagnetic 
ratio, (j> is a phase constant, and 

n(B') = {5[B'-B(r)}), (4) 

is the probabilty of finding a local magnetic field B in 
the z-direction at a position r in the x-y plane. 

Application to studies of the vortex state 

The internal magnetic field distribution n{B) mea- 
sured by yuSR is similar to what one can measure with 
NMR. However, because the muon has a spin equal to 
1/2, it does not possess an electric quadrupole moment 
and hence is a pure magnetic probe. In contrast, n(B) 
measured by NMR often includes quadrupolar broaden- 
ing and in some cases closely spaced quadrupolar satel- 
lites. At its stopping site, a given muon randomly samples 
n(B). This is because the intervortex spacing is typically 
several orders of magnitude larger than the atomic lat- 
tice spacing. For this reason it is often not necessary to 
know precisely where the muon stops in the sample. 

Although the ^SR signal is recorded in the time do- 
main, a fairly accurate visual illustration of n(B) is pro- 
vided by the Fourier transform of P(t) — often called the 
ySR line shape'. An example of a fiSR line shape for 
YBa 2 Cu3 06.95 is shown in figure^ Included in this fig- 
ure is n(B) of equation §5§ obtained from a fit to the 
asymmetry spectrum, the details of which are described 
later in this article. Despite the Gaussian apodization 
used in generating the Fourier transform, the /xSR line 
shape closely resembles n(B) obtained from the fit to 
the /iSR signal in the time domain. Note that n(B) dif- 
fers somewhat from that expected for an ideal FLL, be- 
cause of FLL disorder and the contribution from nuclear 
dipoles. 

As shown in figure [3J the minimum size of the sample 
required for /zSR studies has greatly decreased over the 
years. Today's capability to probe smaller samples comes 
from several key developments. Like neutron scattering, 
the quality of /iSR spectra are limited by counting times 
and background contributions. Shorter data collection 
times and smaller background signals have been achieved 
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FIG. 2: Fourier transform of the TF-/iSR signal from sin- 
gle crystal YBa 2 Cu 3 06.95 at T = 2.5 K and H = 4.94 kOe 
(black curve) computed using a Gaussian apodization func- 
tion with a width of 3 /xs -1 . The peak at B « 4.94 kG is a 
background signal coming from muons stopping outside the 
sample. Also shown is the internal magnetic held distribution 
n(B) (red curve) obtained from fits to the (iSR precession sig- 
nal. The fit assumes n(B) of an ideal FLL (blue curve) gen- 
erated from a Ginzburg-Landau model described later, and 
additional broadening by FLL disorder and nuclear dipole 
moments. 
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FIG. 3: Evolution of the minimum sample size required for 
/iSR studies of the vortex state. 



through the availibility of higher luminosity muon beam 
lines, which allow more of the beam to be focused on the 
sample. Major advancements in /iSR spectrometer de- 
sign [ll|, [13] have also greatly helped in suppressing the 
background signal arising fom muons that miss the sam- 
ple and stop elsewhere. In many setups this is achieved 
by placement of the muon and positron counters within 
the same cryogenic environment as the sample. This 
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tighter geometry makes it easier to restrict the good de- 
cay events to those associated with muons stopping in 
the sample. The background peak that contributes to 
the /LiSR line shape in figure [2] constitutes ~ 12 % of 
the total signal, and is accounted for by adding a polar- 
ization function exp(— a 2 tg t 2 /2) cos^-Bbg* + </>b g ) with 
12 % amplitude to equation ([3]). This /iSR line shape 
was recorded in 1992. As shown later, today this same 
measurement can be done with virtually no background 
signal. However, it is often helpful in the anaylsis of such 
data to have some background, as it provides an accu- 
rate determination of the average magnetic field at the 
sample position. 



THE MUON DEPOLARIZATION RATE 

In the London model the local magnetic field for a 
perfect FLL is 



B(t) = B 
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where the sum is over reciprocal lattice vectors G of the 
FLL. This equation is considered valid provided A>( 
and B <C B C 2 , where £ is the superconducting coherence 
length and B c2 is the upper critical field. The second 
moment of the internal magnetic field distribution n(B) 
corresponding to equation (JSJ) is [13| 



((AB) 2 ) = ((B(v)~(B)) 2 )= V $ 2 X- 



(G) 



where (B) is the average of n(B), <E>o = hc/2e — 2.07 x 
1CT 15 T-m 2 , and r\ = 0.00371 and r\ = 0.003868 for a 
hexagonal and square FLL, respectively [14| . In the Lon- 
don model ((AB) 2 ) is independent of the applied field. 
In high-T c superconductors, where the effective mass of 
the charge carriers is larger in the direction perpendicu- 
lar to the CuC>2 layers compared to that parallel to the 
CuC>2 layers, A is highly anisotropic. If the field is ap- 
plied along the c-axis, then A = (AaA/,) 1 / 2 = X a b, where 
X a b is the in-plane magnetic penetration depth describ- 
ing screening currents flowing in the d-b plane around a 
vortex. Applying the field perpendicular to the c-axis re- 
sults in the formation of Josephson vortices sandwiched 
between Cu02 layers. In this case the screening currents 
flow in the a-c or b-c plane, and A = (AqAc) 1 / 2 = X ac or 
A = (AfcAc) 1 / 2 = Abe, respectively. The out-of-plane pene- 
tration depth A c is typically much greater than 5000 A, 
so that ((AB) 2 ) is significantly reduced in this geome- 
try. For such large values of the magnetic penetration 
depth, the width of the /j,SR line shape in the vortex 
state becomes comparable to that in the normal state, 
and this prevents an accurate determination of the tem- 
perature dependence of A. The absolute values of A a , 
Xb and A c cannot be isolated by /iSR without a priori 
knowledge of the mass anisotropies. For example, if it is 



known that Ah = A Q , then A a can be determined from a 
measurement of A a f>, and A c can be determined from a 
subsequent measurement of A ac . This has been done for 
the heavy-fermion superconductor UPt3 jl5j |. 

For polycrystalline samples an often used approxima- 
tion is to assume n(B) is a simple Gaussian distribution 
of fields of width a 



n(B) 
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in which case equation ([3]) has the following form 

P x (t) = exp(-o- 2 t 2 /2)cos(-f^(B)t + 0). (8) 

The width a is thus equivalent to the damping rate of 
the asymmetry spectrum, called the 'muon depolariza- 
tion rate'. There is generally a temperature-independent 
contribution to a from randomly oriented nuclear dipole 
fields, which can easily be determined from measure- 
ments above T c . The resultant measured value of a can 
then be used to obtain the second moment of the Gaus- 
sian distribution n(B) 



((AB) 2 ) = a 2 hi 



(9) 



By equating equation © to equation ©, one gets an 
estimate of A 



(10) 



For large anisotropy 7 = [ra c jm a b) x ^ 2 = X c /X a b > 5 the 
effective penetration depth A is primarily determined by 
X a b, such that 



X — fXab ■ 



(11) 



where / = 1.23 [16] . For 7 = 1 to 5, / varies from 1 to 
1.23. 

Early /iSR measurements of a were used to deduce 
information about 1/A 2 & using equations (fit))) and (fTTjl . 
From observations of a weak temperature dependence for 
a at low T it was concluded that the pairing state symme- 
try of high-Tc superconductors is s-wave [IT], E, 0, [2(| • 
A linear scaling of T c with a was also found (the so-called 
'Uemura plot'), which has been interpreted via equation 
(JTU]) to imply T c oc 1/A 2 b , which in turn is proportional to 



the superfluid density n s [21|, |22| . While the dependences 



of a for newly discovered superconductors continue to be 
reported in the literature, the assumed relation ctoc 1/A 2 fc 
is not precisely correct. One reason is that the use of 
equation (fTT))) assumes one is aware of and understands 
all contributions to the measured value of a. Brandt [l3| 
has shown that for straight rigid vortex lines, a is in- 
creased by both pinning-induced random displacements 
of the vortices from their ideal position in the FLL and 
thermal fluctuations of the vortex lines. On the other 
hand, the layered nature of the high-T c superconductors 



5 



means that the vortices can be highly flexible. In this 
case segments of a flux line are susceptible to pinning or 
thermal-induced displacements, which reduce a. In ei- 
ther situation the size of the effect will be dependent on 
temperature and the applied magnetic field. There may 
also be contributions to a from static electronic moments. 
This is likely to be the case in rare-earth and lightly- 
doped cuprate superconductors. Even if identified, these 
contributions to a cannot be reliably separated from the 
depolarization rate associated with A a fc. 

A second consideration is that Yaouanc, Dalmas de 
Reotier and Brandt [23| have shown that the second mo- 
ment due to an ideal hexagonal FLL is more accurately 
given by 



((AS) 2 ) = 0.00371$ 2 A- 4 /^), 



(12) 



where b = B jB c i is the reduced field and f v (b) is a uni- 
versal function that accounts for the finite size of the 
vortex cores. It is often assumed that the vortex cores of 
high-T c superconductors are small (i.e. « 15 A), so 
that little error is introduced in neglecting them. How- 
ever, measurements by other techniques indicate that the 
vortex core radius can greatly exceed 15 A in under- 
doped and overdoped samples 24], 25, 26}. As explained 
in Refs. [23,|27|, f v is strongly dependent on b. Moreover, 



Brandt [27] has shown that the often used approximation 
f v = 1 is really only satisfactory over a narrow field range 
and only for k > 70. The latter condition is not always 
satisfied in high-T c superconductors. 

A third problem is the assumption that n(B) is a Gaus- 
sian distribution. In a polycrystalline sample n(B) is an 
average over all orientations 13, 28}. An example of what 
n(B) should look like in an unoriented powder sample 
is shown in figure 0] The fact that many of the early 
/iSR line shapes from measurements on polycrystalline 
samples of high-T c superconductors looked Gaussian-like, 
suggests that the samples were of poor quality with a 
high-degree of FLL disorder. 

Given these issues, the dependences of a should be 
viewed as providing the correct qualitative trends of 1/A 2 
only. Indeed the Uemura plot was an important achieve- 
ment, as it established for the first time a qualitative 
relationship between T c and the superfluid density. How- 
ever, the current state of affairs regarding the search for 
a microscopic theory of high-T c superconductivity is such 
that precise quantitative relations are desired. In recent 
years, experimental studies of 1/A 2 have shown that T c 
is a sublinear function of n s in the underdoped regime 
29l . 3oL 3l| and that there is no single universal Uemura 
relation for high-T c cuprates [32I ]. Likewise, experiments 
on single crystals (33|, including fiSR [34, 35 1, have estab- 
lished a limiting low-temperature linear T dependence of 
X a b, that is not evident in the earlier measurements of 
a. This behaviour is consistent with a d x 2_ y 2-w&ve su- 
perconducting order parameter, which has been clearly 
identified by phase-sensitive techniques [36] ]. 
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FIG. 4: Theoretical field distributions for an unoriented poly- 
crystalline sample calculated from the 'Kossler model' for the 
powder average of an anisotropic superconductor [2|| and in- 
cluding FLL disorder. The parameters used to generate the 
blue curve come from the fit of single crystal YBa2Cu306.95 
in figure [2] The red curve is calculated with a higher de- 
gree of FLL disorder (i.e. ad of equation (|15[) is three times 
larger), and is closer to a Gaussian distribution. (Note that 
the two curves for n(B) are not normalized with respect to 
each other.) 



FULL LINE SHAPE ANALYSIS 

To extract values of A and the vortex core size from 
TF-/1SR measurements, one must assume a model for the 
internal magnetic field distribution n(B). This in fact is 
the largest source of uncertainty in these kinds of mea- 
surements. To date, the analysis of /iSR measurements in 
the vortex state have primarily relied on models based on 
either the phcnomenological London or Ginzburg-Landau 
(GL) theories. The reason is that in some limiting cases, 
analytical models for the spatial field profile B(r) exist 
for these theories. While the analytical London and GL 
models provide very good fits to ^SR measurements of 
the FLL in real materials, neither are generally valid, 
and considerable effort has been devoted in recent years 
to understanding the meaning of the fitted parameters 
A and £. In particular, it is important to consider the 
following cautionary statements, which are discussed in 
further detail in the following sections: 

• Generally speaking, the fitted parameter A is 
not a measure of the magnetic penetration 
depth. This is because A is influenced by devi- 
ations of the assumed theoretical model for n(B) 
from the actual field distribution of the FLL in the 
sample. This point was made in Refs. 23, 3?], S, 
H, O, EH, S3, EH, 0, EU . For this reason, the mag- 



netic penetration depth actually corresponds to the 
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extrapolated H^O value of A. 

• Likewise, the fitted parameter £ is generally 
not a measure of the superconducting coher- 
ence length. This point is explained in detail in 
Ref. [i^l ■ Instead £ qualitatively tracks changes in 
the vortex core size that have nothing to do with 
the size of the Cooper pairs. 

On the theoretical front, great strides have been made 
in numerically calculating n{B) from the microscopic 
equations. The predictions of the microscopic theory 
provide one means of determining the accuracy of the 
phenomenological approach to modelling fiSK measure- 
ments. 



Finite core size 

Modified London models. 

The London model does not account for the spatial 
dependence of the superconducting order parameter and 
consequently equation §5§ breaks down at distances on 
the order of £ from the vortex core centre — i.e. B(r) 
diverges as r — > 0. To correct this, the G sum can be 
truncated by multiplying each term in equation ([5]) by a 
cutoff function F(G). It is important to stress that the 
appropriate form of F(G) depends on the precise spatial 
dependence of the order parameter in the the vortex core 
region, and this in general depends on temperature and 
magnetic field. 

Analytical expressions for F(G) exist only for certain 
limiting cases. Some time ago Brandt suggested the use 
of a smooth Gaussian cutoff factor 



J F 1 (G)=exp(-aG 2 ^ 2 ). 



(13) 



If there is no dependence of the superconducting coher- 
ence length £ on temperature and magnetic field, then 
changes in the spatial dependence of the order parame- 
ter around a vortex correspond to changes in a. By solv- 
ing the GL equations, Brandt determined that a = 1/2 at 
fields near B c2 |X3|, [46|, l47|] , and arbitarily determined that 
a w 2 at fields immediately above B c \ [HI . More accu- 
rate values of a have since been obtained from precision 
solutions of the GL equations. For an isolated vortex in 
an isotropic extreme (k^>1) s-wave superconductor, a is 
found to decrease smoothly from a = l at B c \ to a«0.2 
at B C 2 [1^1 . Even so, it is rather presumptuous to in- 
corporate the field dependence of a into the analysis of 
/iSR spectra, as the vortices in real materials are not iso- 
lated. Furthermore, the additional effects of temperature 
on the spatial dependence of the order parameter in the 
vortex core region are not accounted for in calculations 
of a(B). For example, Laiho et al. 4^, H3| have shown 
by comparison to solutions of the quasiclassical Eilen- 
berger equations for a d x 2_ y 2-w&ve superconductor that 



a is temperature dependent. More recently they have de- 
termined the magnetic and temperature dependences of 
the cutoff function used in the London model by numer- 
ically solving the quasiclassical Eilenberger e quat ions for 
the vortex state of an s-wave superconductor [44| . These 
calculations show that the London model with the proper 
cutoff function provides a reasonable description of the 
field distribution of the FLL in real superconductors. 

Generally, a is fixed when using the modified London 
model to fit /iSR spectra, so that the temperature and 
magnetic field dependences of the order parameter are 
indicated by variations in the fitted value of £. Such a 
procedure has been used to account for the finite size of 
the vortex cores in a number of /iSH studies of the vortex 
state in type-II superconductors 0, HI, 52, SHU!. 
Unfortunately, none of these studies were done at suffi- 
ciently low enough temperature to make a direct com- 
parison of the field dependence of £ with the theoretical 
predictions for a(B) at T — K. 



Gmzburg-Landau models. 

In recent years, modified London models for B(r) have 
been abandoned by some in favour of models based on 
GL theory. The appealing aspect of the GL models is 
that the spatial variation of the order parameter is nat- 
urally built into the theory. The drawback is that GL 
theory assumes that the order parameter varies slowly in 
space and is strictly valid only near T c . Despite these 
limitations, GL theory has proven to be highly success- 
ful in describing variations of n(B) as measured by /iSR, 
yielding accurate quantitative values of A and £ in certain 
cases. As is the case in using modified London models, 
the key is to be careful with the interpretation of the 
fitted values. 

The GL equations for the ideal Abrikosov vortex lat- 
tice can be solved by a variational method [56j]. At low 
reduced fields b — Bj B C 2 <C 1 and for k>1, an excellent 
analytical approximation to the spatial field profile in GL 
theory is [23j 

B{v)=B a {l-b±)Y j e lG ^ (G) ' (W) 



G 



where F(G)=u K x (u) , u 2 = 2£ 2 G 2 (1 + b 4 ) [1 - 25(1 - bf] , 
and K\ (u) is a modified Bessel function. Note the cutoff 
function F(G) depends on the local internal magnetic 
field B. 

Brandt later introduced an interative procedure for 
solving the GL e quat ions for arbitrary 6, k and vortex- 
lattice symmetry [57[. While more generally applicable, 
fitting fiSK spectra with this iterative method requires 
a tremendous amount of computer time. Fortunately, it 
does not seem that this is necessary. In a recent ^iSR 
study of the 1ow-k type-II superconductor V, qualita- 
tively similar results were obtained from analyses with 
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the variational and iterative solutions of the GL equa- 
tions 



42] 



Vortex symmetry. 

In the above models the individiual vortices are as- 
sumed to have circular symmetry. However, this is not 
the case in real type-II superconductors with anisotropic 
Fermi surfaces and/or superconducting energy gaps 58|, 



60, 61 



|62j. At magnetic fields where there is an 
appreciable degree of overlap of the vortices, the FLL 
adopts the symmetry of the individual flux lines. In cases 
where the deviation from circular symmetry is extreme, 
the symmetry of the FLL can be accounted for in the 
model for B(r). This has been successfully done in /xSR 
studies of YNi 2 B 2 C [H^L LuNi 2 B 2 C [H, Sr 2 Ru0 4 0, 
V 3 Si [13| and Nb 3 Sri ij65|. In the case of V 3 Si, it was 
possible to extract from the /iSR measurements a grad- 
ual hexagonal-to-square transformation of the FLL sym- 
metry in good agreement with STS imaging experiments. 
A couple of measurements from this study are shown in 
figure [3 

In Refs. (HiL [55L the theoretical models used for 
analysis of the TF-//SR spectra are based on London 
theory [(3(J [6l|. London models that account for the 
fourfold symmetry of the ci^^^-wave order parameter 
in cuprate superconductors have also been developed 
[371 16a ] . Anisotropy may also be incorporated into ana- 
lytical GL models. The TF-^SR study of Sr 2 Ru0 4 (si 
employed a GL model that accounts for p-wave symme- 
try and Fermi surface anisotropy (67l . [6a ]. Since all of 
the theoretical models accounting for anisotropy intro- 
duce additional fitting parameters, they are seldom used 
in /iSR studies. When anisotropy is not included in the 
model for B(r), the fitted values of A and £ are angle- 
averaged length scales. 



Vortex Lattice Disorder 

3D vortex lines. 

The effects of random pinning and thermal fluctuations 
of the vortices on n{B) depend very much on the rigidity 
of the vortex lines. In general, reliable information on the 
length scales A and £ is obtained only when the vortices 
are fairly straight and parallel. In this case, pinning or 
fluctuation induced displacements of the vortices from 
their positions in the perfect FLL increase the width of 
n(B). The effect is nearly equivalent to smearing n(B) 
by convolution with a Gaussian distribution of fields 47 j 



n(B) = 



1 



2-nod 



exp 



B-B' 



n{B')dB' , 

(15) 



where ad is related to the root mean square displacement 
(wf) 1 / 2 of a vortex line about its average position. From 
/iSR studies on numerous superconductors with 3D-likc 
vortices it has been determined that ad cx B/X 2 ab . This 
simply means that stronger overlap of neighbouring vor- 
tices reduces the degree of disorder. 

The above treatment is a reasonable approximation 
of the weak FLL disorder in the 3D 'Bragg-glass' phase 
of a type-II superconductor, in which quasi-long-range 
translational order and perfect topological order are pre- 
served |7Qj . Note that a more accurate expression for the 
variance of the FLL field distribution in the Bragg-glass 
phase has been derived 7l|, but has yet to be used for 
the analysis of /xSR measurements. More recently, it has 
been shown that when disorder is strong enough to pro- 
duce a 3D 'vortex-glass' phase, in which the topological 
order of the FLL is not retained, the ySR line shape is 
slightly skewed in the opposite way with a 'tail' on the 
low-field side 72, 73j. In this situation, the disorder can- 
not be handled in a trivial way. 



2D pancake vortices. 

The vortex lines in layered superconductors are gener- 
ally considered to be comprised of coupled 2D 'pancake' 
vortices. When the interlayer coupling is weak, the vor- 
tex lines are very soft, and pinning and/or thermal fluc- 
tuations of the pancake vortices shorten the high-field 
tail of n(B) 0[zilz|- 

There has been some debate on how disorder of the 
FLL should be treated in /iSR studies of high-T c su- 
perconductors. An early /iSR study by Harshman et 
al. [7H ] arrived at the conclusion that point distortions 
of the vortex lines was significant in highly anisotropic 
Bi 2 Sr 2 CaCu 2 0§+5, as evidenced by narrow symmetric 
/iSR line shapes at low temperatures and high mag- 
netic fields. On the other hand, under similar condi- 
tions broad asymmetric /LtSR line shapes were observed 
for YBa 2 Cu30j / , consistent with vortex lines that do 



not wander significantly along their length. The dimen- 
sionality of the vortex lines is dependent on the ratio 
Js/Xab, where 7 is the mass anisotropy and s is the 
spacing between Cu0 2 planes. The observations re- 
ported in Ref. [lH are consistent with the fact that 7 
is nearly two orders of magnitude smaller in optimally 
doped YBa 2 Cu30 !/ than in Bi 2 Sr 2 CaCu 2 0s+5. Thus in 
high-r c cuprates the picture of weakly coupled 2D pan- 
cake vortices should be viewed as a limiting case. 

Recently, Harshman et al. [77| have asserted that point 
distortions of the vortex lines in ultra-clean samples of 
YBa 2 Cu 3 O y are as important as in Bi 2 Sr 2 CaCu 2 08+5- 
By fitting the total second moment of the /j,SR line 
shape for fully doped YBa 2 Cu 3 07 to a phenomenolog- 
ical model containing two independent parameters for 
pinning-induced distortions of the FLL, Harshman et al. 
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FIG. 5: Fourier transforms of the TF-/iSR signal from single crystal VsSi at T = 3.8 K for applied magnetic fields of H = 30 kOc 
and H = 50 kOe directed parallel to the [001] axis (from data in Ref. 64]). The solid red curve in each panel is a Fourier 
transform of the fit in the time domain (Note: because of apodization, this is slightly different than the fitted n(B), an example 
of which was shown in figure [2}. Also shown are contour plots of the field profile B(r) obtained from the fits to a GL model. 
The results are in good agreement with STS measurements on VsSi [691 ] . 



argued that X a b(T,H) is consistent with s-wave super- 
conductivity. However, there are several problems with 
the analysis of the /iSR data in Ref. [77| ■ Most notably, 
the second moment of n(B) is assumed to be given by 
equation ([6]), which does not account for the finite size 
of the vortex cores. Also, both the determined value 
« = A ah /Ca6 = 43.8± 1.8 and the lowest field (5 = 0.05 T) 
considered in Ref. 77| are too small for equation ^ to 
apply. The smaller value of the second moment found 
at B = 0.05 T in Fig. 2 of Ref. [77] is in fact expected, 
because at this field ((AS) 2 ) is more appropriately de- 
scribed by equation (12) of Ref. [27J 



hK 2 * 2 

((ABO - 



(16) 



in YBa2Cu30 y and 
more rigid than 



Although the vortex lines 
La2_ a; Sr a ;Cu04 are certainly more rigid tnan in 
Bi2Sr2CaCu208+5, they do soften in underdoped sam- 
ples due to an increase in anisotropy 7. Experiments 
on oxygen-deficient YBa2Cu30 y thin films O in- 
dicate that for applied magnetic fields currently attain- 
able in a TF-/iSR experiment (i? = 80 kOe), the FLL 
remains 3D at low T, except perhaps in lightly doped 
samples. However, the dimensionality of the vortex sys- 
tem depends on the degree of sample disorder, which can 
cause pinning-induccd misalignment of the pancake vor- 
tices that make up a vortex line. Since YB^CuaOy thin 
films can be strongly disordered, these experiments are 
sample dependent and don't necessarily tell the story in 
clean samples. In fact, recent mutual inductance mea- 
surements on severely underdoped thin films indicate 
that YBa2Cu30j, is really quasi-2D only near T c [30l |. 



The weak field dependence of the Josephson plasma res- 
onance in YBa2Cu306.5o single crystals at low T is fur- 
ther evidence that the vortices maintain a 3D charac- 
ter in clean underdoped samples [80l ]. As for under- 
doped La2- a; Sr a: Cu04, it has been argued from a com- 
bined /iSR and small-angle neutron scattering study of 
Lai gSr 1CUO4 [72| | that the vortex lines remain fairly 
rigid at low T. 



VORTEX CORE SIZE 

Superconductivity is strongly suppressed in the vortex 
core. A loose definition of the vortex core size is that it 
is a region of radius r where £ is the characteristic 
length scale for spatial variations of the superconducting 
order parameter ip — i.e. the GL coherence length. How- 
ever, this definition is not particularly satisfactory, espe- 
cially since GL breaks down at low temperatures. Two al- 
ternative definitions of the vortex core size have emerged 
from numerical solutions of the Bogoliubov-de Gcnncs 
and quasiclassical Eilenberger equations. The first comes 
from the initial slope of the pair potential A(r) near the 
centre of the vortex core. Assuming A(r)ocr close to the 
core centre, the core radius may be defined as 



£1 = A / Urn 



A(r) 



(17) 



where Ao ~ l/£o is the bulk superconducting energy gap 
at zero temperature, and £0 is the BCS coherence length. 
As shown in figured, r = £1 is the radius at which the lin- 
ear slope extrapolates to A(r) = Ao- The absolute value 
of the spatial supercurrent density profile \ j(r)\ provides 
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FIG. 6: Spatial variation of (a) the pair potential A(r) and 
(b) the absolute value of the supercurrent density j(r) for an 
isolated vortex, where r is the radial distance from the vortex 
centre. 



a second way of denning the core size. As shown in fig- 
ure [SJa, \j(r)\ exhibits a maximum at a distance r = ro 
from the core centre. 



Magnetic field dependence 

Figure [7J shows the magnetic field dependence of the 
vortex core size determined by /iSR in the k>1 super- 
conductors VsSi, NbSe2 and LUM2B2C. Similar results 
have been reported for CeRu2 [53|, YNi 2 B 2 C [54( and 
NbsSn 65]. In agreement with calculations from the qua- 
siclassical Eilenberger equations 81 , [SH , r and £ a b ex- 
hibit qualitatively similar dependences on magnetic field. 
As explained in Ref. the magnetic field dependence 
of ro is partially due to the 'vortex squeezing effect' — 
which refers to the increasing overlap of the j(r) profiles 
of individual vortex lines with increasing H (see figure 2 
of Ref. O). 

It is important to stress that the field dependence of 
the parameter £ a b that comes from fits of /xSR measure- 
ments using the models described in section 4.1, does 
not mean that the superconducting coherence length is 



changing with field. Rather, it simply means that the 
cutoff factor in these models significantly varies due to 
changes in the slope of A(r) in the vortex core region. 
For this reason the behaviour of £ a ;, is expected to fol- 
low that of £1 of equation (JTTJ), and hence is considered a 
measure of the vortex core size. In general, the maximum 
value of the cutoff parameter £ a & measured by /iSR corre- 
sponds to the GL coherence length calculated from H C 2 
(i.e. £gl = (Qo/ZTrH^) 1 / 2 ). At low fields, where the vor- 
tices are weakly interacting, the fitted value of £ a b more 
or less agrees with that expected from H C 2 (see figure [7]). 

The observed shrinkage of the vortex cores at higher 
magnetic fields qualitatively agrees with theoretical cal- 
culations for high-K superconductors done in the frame- 
work of the quasiclassical Usadel [H 0] and Eilenberger 
equations [HI, |s3|- These theoretical works show that 
changes in the spatial variation of j(r) (and hence B(r)) 
and A(r) with increasing magnetic field are directly re- 
lated to changes in the electronic structure of the vor- 
tex cores that occur with increased vortex-vortex inter- 
actions. Physically, the situation is analogous to the 
electronic band structure that forms when atomic or- 
bitals are combined in a crystal lattice. When atoms 
come together to form a lattice, the wave functions of 
the electrons anchored to individual atoms overlap the 
wave functions of the electrons of neighbouring atoms. 
This results in the formation of energy bands, which in 
the case of a metal, leads to delocalized electrons that are 
capable of carrying electric and thermal currents. In a 
similar way, as vortices are brought closer together by in- 
creasing the applied magnetic field, there is an increased 
overlap of the wave functions of the single-particle ex- 
citations (quasiparticles) bound to an individual vortex 
[83l | with the quasiparticle states of neighbouring vor- 
tices [8ll [82], [3 Hlf , resulting in the intervortex transfer 
of quasiparticles. Since the higher-energy bound quasi- 
particle core states are more spatially extended, they 
delocalize first. The reduction of the vortex core size 
with increasing magnetic field is thus understood as be- 
ing due to the sequential delocalization of quasiparticle 
core states, starting with the highest-energy bound state, 
and terminating with the complete delocalization of the 
lowest-energy bound state. 

While this picture is strongly supported by a compar- 
ison of /iSR measurements of the vortex core size with 
the data from other techniques that probe electronic ex- 
citations, Kogan and Zhelezina (8?| have proposed an al- 
ternative theoretical explanation for the field dependence 
of the core size that is based on weak-coupling BCS the- 
ory. Their model predicts that in clean high-K supercon- 
ductors, the intrinsic superconducting coherence length 
is not independent of magnetic field, but actually de- 
creases with increasing field. Recently, DeBeer-Schmitt 
et al. j88j measured an unusual field-independent FLL 
form factor in CeCoIns by small-angle neutron scatter- 
ing, suggestive of a field-dependent coherence length (see 
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figure [5]). Even so, there may be another explanation for 
their data that warrants further experimental investiga- 
tion of the FLL form factor in other materials. A seem- 
ingly direct test of the Kogan-Zhelezina theory would be 
to study the magnetic field dependence of the core size 
in impurity-doped superconductors. Unfortunately, the 
addition of impurities is likely to interfere with the inter- 
vortex transfer of quasiparticles and add disorder to the 
FLL. Thus the interpretation of such experiments would 
probably be ambiguous. 

Specific heat 

Specific heat measurements are sensitive to both lo- 
calized and delocalized quasiparticle excitations. In the 
vortex state the specific heat is usually described by 

c(T, B) = 7 (T, B)T + aT 2 + f3T 3 , (18) 

where the last term /3T 3 is the phonon contribution. The 
first term is the electronic contribution to the zero-energy 
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quasiparticle density of states (DOS) N(E = 0), i.e. the 
spatial average of the local DOS N(E = Q, r) at the Fermi 
level E = 0. The 'Sommerfeld coefficient' is given as 



7(£) = lim T ^ c(T,B)/TocN(0) . 



(19) 



The second term in equation (jT5J) is not always clearly 
present in experimental data, and only recently has the- 
oretical progress in understanding its origin been made. 
Nakai et al 89] have shown there are two electronic con- 
tributions to the aT 2 term. The first comes from low- 
energy states (E > 0) that appear near the Fermi en- 
ergy at B 7^ 0. The second contribution comes from the 
'Kramer-Pesch effect', which is a shrinking of the vor- 
tex cores with decreasing temperature that is associated 
with the thermal depopulation of the bound quasiparticle 
core states. As discussed elsewhere 0, Hoj], the Kramer- 
Pesch effect has been measured by fiSH in a variety of 
superconductors. 

In an isotropic-gapped superconductor, the low-lying 
excitations a re g enerally considered to be confined to the 
vortex cores [83(. In this case j(B) is proportional to the 
product of the area of a vortex core 7r£ 2 and the density 
of vortices 90], so that 



7(B)CX7T^ 2 B. 



(20) 



Thus, if the size of the vortex cores is independent of mag- 
netic field, then j(B)(xB. In highly anisotropic-gapped 
superconductors, including those with line or point nodes 
at the Fermi surface, the low-lying excitations near the 
gap minima become the dominate contribution to j(B). 
Volovik showed that for the case of a <i-wave superconduc- 
tor, the superfluid flow around a vortex lowers the energy 
of quasiparticles delocalized near the gap nodes, leading 
to a finite DOS at the Fermi energy and a corresponding 
nonlinear field dependence for 7(B) [9l|. For the case of 
a highly anisotropic energy gap without nodes, such be- 
havior is expected only if the energy shift exceeds the gap 
minimum. Until recently, this is more or less how spe- 
cific heat measurements on type-II superconductors as a 
function of magnetic field were interpreted. j(B) cx B be- 
haviour was indicative of conventional s-wave supercon- 
ductivity, and a sublinear dependence of 7(B) on field, 
usually 7CX \[B, meant unconventional superconductiv- 
ity. 

In 1999 it became clear that the delocalization of quasi- 
particles brought about by vortex-vortex interactions 
shows up in the field dependence of the specific heat of 
s-wave superconductors — often mimicking the behaviour 
expected for an unconventional pairing-state symmetry. 
Theoretically it was shown from solutions of the qua- 
siclassical Eilenberger equations at T — 0.5T C that the 
zero-energy DOS per vortex N(Q)/B depends on mag- 
netic field [81]. In the same work this unexpected be- 
haviour was shown to be due to a shrinking of the vortex 
cores with increasing field, associated with an increased 



overlap of the zero-energy DOS of one vortex with that 
of its neighbours. This idea was experimentally verified 
in the same year by relating the field dependence of the 
vortex core size measured in NbSe2 by /iSR [H to the 
field dependence of the specific heat of this material [92| • 
More recent calculations in the framework of the qua- 
siclassical Eilenberger theory [93[ show that for an s- 
wave superconductor at T = 0.1T C , the isolated vortex 
behaviour 7(B) cx B is in fact realized at low tempera- 
tures up to a crossover field B*. Above B* , ~f(B) has a 
sublinear dependence on field due to the overlap of the 
low-energy quasiparticle core states from neighbouring 
vortices. Even so, it is rare to find superconducting ma- 
terials that exhibit 7(B) <xB behaviour at low magnetic 
fields, because any anisotropy of the superconducting en- 
ergy gap or Fermi surface reduces the value of B* . This 
also means that the vortex core size measured by /zSR 
will most often be dependent on field. Recently, this was 
shown to be the case even at low field in the marginal 
type-II superconductor V 0]. The field dependence of 
the core size in pure V cannot be explained by the theory 
of Ref. [87|, which does not apply to low- re superconduc- 
tors. On the other hand, an explanation in terms of 
delocalized core states is compatible with the sublinear 
dependence of j(B) on B found immediately above H c \ 
in the low-Ac type-II superconductor Nb [9^ ]. 



Thermal conductivity 

In contrast to the electronic contribution to the spe- 
cific heat in the mixed state, the electronic contribution 
to heat transport comes entirely from extended or delo- 
calized quasiparticles. In superconductors with extreme 
gap anisotropy, such as LuNi 2 B 2 C [HI or YBa2Cu30 y 
[96(, the dominant contribution to the field dependence 
of the thermal conductivity at low T is the Doppler shift 
of the quasiparticle spectrum outside the vortex cores due 
to the circulating supercurrents 97 , 98| . However, for an 
s-wave superconductor, the field dependence of the elec- 
tronic thermal conductivity K e obtained by extrapolating 
measurements to T — > K is a direct measure of the de- 
localization of quasiparticle states bound to the vortex 
cores. 

In recent years, /iSR measurements on s-wave super- 
conductors have established a direct correlation between 
the field dependence of the vortex core size and the 
delocalization of quasiparticle core states as measured by 
thermal conductivity. By imposing the following simple 
model, excellent agreement between these two very 
different kinds of measurements has been demonstrated 



The electronic thermal conductivity n e is proportional 
to the zero-energy density of delocalized states A^jeioc(O), 
and hence grows as the zero-energy density of bound 



12 



cores states N(0) decreases. Since N(0) is proportional 
to tt£ 2 B, then 

K e cx N dcloc cx K(B cl ) 2 - tt((B) 2 )B , (21) 

where tt£,(B c i) 2 and 7r£(i? c i) 2 are the areas of the vortex 
core at B = B cl and B > B cl , respectively. 



V3S1: a good standard. 

In 2004, a /zSR study of n(B) in the mixed state of 
the high-K s-wave superconductor VsSi was reported 64 1 . 
The results of this experiment provide a good standard 
for comparison of /iSR measurements on more compli- 
cated systems. The reason is that at low field, the elec- 
tronic states remain fairly well localized within the vortex 
cores. Consequently, K e (T^0, H) increases slowly with 
increasing magnetic field [9^], c(T — >0, H) is dominated 
by localized electronic states at low H, and the vortex 
core size determined by /iSR is essentially independent 
of H up to H sa 7.5 kOe. As shown in figure [SJi, equa- 
tion (f2"Tj) accurately describes the thermal conductivity 
data using the /xSR values of the vortex core size £ a &. 



V and Nb: marginal type-II superconductors. 

As mentioned earlier, the vortex core size in the low- 
k superconductor V was recently measured by /^SR [42j ]. 
With the magnetic field applied parallel to the (111) di- 
rection of a V single crystal, the vortex core size was 
found to decrease immediately above the lower critical 
field H c \ « 1 kOe. The high degree of derealization of 
quasiparticle core states necessary to cause the observed 
shrinking of the vortex core size at such low field is a con- 
sequence of the large superconducting coherence length 
(~ 300 A) and 10 % to 20 % gap anisotropy. 

While there are no reported measurements of K e (T — > 
0, H) for V , the re are for the related elemental supercon- 
ductor Nb [l00| . T he ele ctronic thermal conductivity of 
Nb reported in Ref. |l00j behaves much like that of VaSi, 
exhibiting a weak exponential increase above H c \. This 
is clearly at odds with the /iSR results for V. However, 
it is also at odds with the specific heat measurements 
of Nb reported in Ref. [94]], which indicate an appre- 
ciable overlap of the quasiparticle core states of neigh- 
bouring vortices immediately above H c i . As argued in 



Ref. [94|, this discrepancy is due to the way in which 



these experiments were done. When the electronic spe- 
cific heat of Nb was measured under zero-field cooled con- 
ditions, no appreciable contribution from localized or de- 
localized quasiparticle states was observed at fields near 
H c \. This i s con sistent with the thermal conductivity 
data of Ref. llOOj measured in monotonically increasing 
and decreasing magnetic field. On the other hand, un- 
der the same field-cooled conditions implemented in the 



/iSR study of V, the specific heat measurements are con- 
sistent with delocalized quasiparticle states immediately 
above H c \. These experiments show that the intervor- 
tex transfer of quasiparticles is disrupted when the FLL 
is highly disordered. While a well-ordered FLL forms in 
the sample under field-cooled conditions, flux entry into 
the sample is impeded by the Bean-Livingston barrier 
under zero-field cooled or monotonically increasing field 
conditions. The field dependence of n e (T^0, H) has yet 
to be measured in Nb or V under field-cooled conditions. 



NbSe2: multi-band superconductivity. 

For a long time the mixed state of NbSe2 was con- 
sidered an archetype of the FLL in conventional type-II 
superconductors. However, over the past 6 years a num- 
ber of experiments on NbSe2 have provided strong evi- 
den ce for dis tinct energy gaps on different Fermi sheets 
H [lOl EqI • While there is some debate over whether 
the sm aller energy gap resides on the Se or Nb Fermi 
sheets 1031 ] . the experimental signatures of multi-band 
superconductivity are similar to those of MgB g. 

Extending theoretical work on MgB 2 [lol HoR fiol 
107l ] to the case of NbSe2 , suggests that the nonlinearity 
of j(B) observed at low B can be attributed to the oc- 
currence of a smaller energy gap Ao on one of the Fermi 
surfaces. In these theories the electronic states of a vor- 
tex are dependent on the combined contributions from 
two different bands, characterized by large and small su- 
perconducting energy gaps and different Fermi surface 
anisotropics. At low magnetic field the dominant contri- 
bution to the low-energy quasiparticle core states comes 
from the small-gap band. Because of the small size of the 
energy gap, these core states are loosely bound to the vor- 
tex and easily delocalize via overlap with the low-energy 
core states of neighbouring vortices. Consequently, the 
vortex core is large at low field, and shrinks with increas- 
ing magnetic field. However, in contrast to a single-band 
superconductor, low-energy quasiparticle states continue 
to be localized around the the vortex cores at high fields 
due to the contribution from the large-gap band. In par- 
ticular, the large-gap band results in electronic states 
that are more confined to the vortex core. While the 
electronic vortex core states o f the individual bands can 
be probed separately by STS [l08l ] . this does not imply 
that there are two different kinds of vortices separated in 
real space. Both bands contribute to a single vortex, but 
the magnetic field dependence of each contribution is dif- 
ferent. Moreover, the energy gaps of the two bands are 
related in the sense that only a single superconducting 
transition occurs. 

In a ^SR study of NbSe 2 at T = 0.02 K, in which ther- 
mal excitations of the bound quasiparticle core states are 
largely frozen out, Callaghan et al. [4l[ explicitly de- 
tected large vortex cores at low magnetic fields (see fig- 
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ure fTO")) . With increasing H the core size rapidly shrinks. 
This behaviour is attributed to the ease at which the 
loosely bound core states delocalize via their overlap with 
the core states of neighbouring vortices. As shown in fig- 
ure HOI the core size saturates above H w 0.1H C 2. At 
these fields the experiment probes the smaller vortices 
associated with the large-gap band. The saturation of 
the core size reflects the reluctance of the tightly bound 
core states of these smaller vortices to delocalize. This 
interpretation of the /iSR measurements is in accord with 
what is known about the vorti ces i n MgB2. At fields far 
below H C 2, Eskildsen 



The unavailability of single crystals of sufficient size has 
so far prevented a study of £, a b(H) in MgB2 by /xSR. 

The two-gap interpretation of the field dependence of 
the core size in NbSe2 is also supported by a remarkable 
agreement with the field dependence of n e (T — > 0, H ) 
measured independently in Ref. [9§] (see figure [SJd) . In 
contrast to VsSi, K e (T — > 0,H) rises rapidly with in- 
creasing field just above H c i, indicating the presence 
of highly delocalized quasiparticle states. This is the 
same behaviour that is observed for the electronic ther- 



et 



MgB2 by STS that have a large core size relative to es- 
timates from H C 2 and do not have localized core states 
This is consistent with the full field dependence of th e 
core size inferred from specific heat measurements llO' 



mal conductivity of MgB2 110] and the heavy- fermion 



108| detected vortices in compound PrOs 4 Sbi 2 (ml lll^, and is understood to 

arise from two or more superconducting energy gaps. As 
is the case for MgB2, intermediate size single crystals 
of PrOs4Sbi2 are not yet available for /iSR studies, and 
hence the field dependence of the vortex core size has not 
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FIG. 10: Magnetic field dependence of the cutoff parameter 
£ a b from fiSK measurements on single crystals of NbSe2 at 
different temperatures [3. l4l|. 
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FIG. 11: The magnetic field dependence of £ a (, in 
YBa2Cu3 06.95 determined by fiSK [381 ] and the normal- 
ized electronic thermal conductivity K e (H) /«; e (0) measured 
in YBa2Cu 3Q6. 9 by Chiao et al. [9(| and in YBa2Cu307 by 
Hill et al. [l!5j . respectively. All data are extrapolations to 
T — > K. Note, the field dependence of shown here is 
a more accurate analysis of the measurements of Ref. [38| . 
In particular, the modified Bessel function K\{u) in equa- 
tion 1)141) was calculated numerically, whereas the analysis 
of Ref. [38| assumed the following analytical approximation: 
R\(u) = (tt /2u) 1/2 exp(-u) , valid for u»l. 



been measured. 



YBa2CusO y : highly anisotropic and multiple gaps? 

As shown in figure [TTJ the vortex-core size in the 
high-T c superconductor YBa2Cu3C>6.95 [38| shrinks with 



increasing magnetic field, such that £ a b(H) — £ a &(0) oc 
H^ 1 / 2 . At first glance this behaviour resembles the pre- 
vious examples, where the shrinking of the vortex cores 
has been associated with the derealization of quasiparti- 
cle core states. However, in YBa2Cu30 y the low-energy 
quasiparticle core states are expected to be extended 
alon g th e nodal directions of the rf x 2_ y 2-wave gap func- 
tion [5^, [82J, and hence are already delocalized at low 
field. Furthermore, the value of £ a b deduced by /zSR 
at low H seems rather large. As discussed in Ref. [ifjj ]. 
the large vortex-core size at low H appears to be unique 
among the cuprates and likely results from proximity- 
induced supercondu ctivit y on the CuO chains. Calcu- 
lations by Atkinson [l 1 3I | show that large two-fold sym- 
metric vortices result from coupling of the CuO chain and 
CuC>2 plane layers (see figure [IS}. The vortex core size 
in YBa2Cu3C>6.95 saturates near H = 40 kOe, at which 
the F LL begins a gradual hexagonal-to-square transition 
114j . The d x 2_ y 2-w&ve superconducting order parameter 



is predicted to produce fourfold symmetry in the spatial 
field profile B(r) around the vortex cores [59(. Thus the 
fourfold symmetry of the FLL at high fields likely results 
from strong overlap of the intrinsic d x 2_ y 2--w&ve vortices. 

A comparison of the /iSR data to electronic thermal 
conductivity measurements by Chiao et al. (9fl | seems 
to imply that there are bound core states which delocal- 
ize with increasing magnetic field (see figure [TTj) . How- 
ever, the behaviour of n e (H) is more likely due to a field- 
induced Doppler shift of the quasiparticle spectrum out- 
side the vortex cores [§3, [98} . More recent measurements 
on u ltraclean single crystals of YBa2Cu307 by Hill et al. 
[l 1 5l | show a saturation of n e (H) at low H, argued to 
indicate an exact cancellation of the contributions from 
the Doppler shift and scattering of the quasiparticles by 
the vortices. Thus, consistent with a highly anisotropic 
gap, there is no established correlation between the size 
of the vortex cores and electronic thermal conductivity 
measurements on very clean samples. 



RNi2B2C: highly anisotropic or multiple gaps? 

For some time now, experimental evidence has 
mounted that the superconducting energy gap in rare- 
earth nickel borocarbides i?Ni2B2C (R = Dy, Ho, Er, 
Tm, Lu and Y) is highly anisotropic, with line or point 
nodes. For example, thermal conductivity measurements 
on LUNI2B2C reveal the existence of highly delocalized 
quasiparticles above H c i [95| . However, there is also ev- 
idence suggesting that these boroc arbide materials ex - 
hibit multi-band superconductivity |l!6l . 117 . 118 . 119| . 
Figure [9J; shows that unlike YBa2Cu30j,, for LUNI2B2C 
there is good agreement between fj.SH measurements of 
the vortex core size and the field dependence of the elec- 
tronic thermal conductivity. This indicates that there 
are quasiparticle vortex-core states at low temperatures 
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FIG. 12: Contour plot of the spatial dependence of the magnitude of the supercurrent density j(r) for vortices in a singl e-layer 
d x 2_ y 2-wave superconductor (left) calculated from self-consistent solutions of the Bogoliubov de Gennes equations [l 13t ] . Also 
shown are the results of similar calculations for a 2 layer superconductor, consisting of a chain layer coupled to a single d x 2_ y 2- 
wave superconducting plane. The middle panel shows what the vortices look like on the plane layer, and the right panel shows 
what they look like on the chain layer. 



that delocalize with increasing magnetic field. Since the 
anisotropy of the superconducting energy gap manifests 
itself in the electronic structure of the vortices, the su- 
perconducting gap of LuNi2B2C must be less anisotropic 
than in YBa2Cu30j,. Moreover, measurements done to 
date are not inconsistent with multiple superconducting 
gaps. The /xSR measurements on LUM2B2C by Price 



et at. [55j only extend down to T = 2.5 K. As in the 



case of NbSe2, a clear signature of multiple-gap super- 
conductivity might be seen in the field dependence of 
the vortex-core size at lower temperatures, where ther- 
mal excitations of the quasiparticle core states are largely 
frozen out. Interestingly, /iSR measurements of S, a b{H) 
in YNi 2 B 3 C at T = 3 K by Ohishi et al. f54j display a 
clear saturation above H ss 5 kOe indicative of a second 
superconducting gap. 



Magnetic field dependence of \ a b 

It is important to stress that the field dependence of 
the fit parameter A a & generally does not imply that the 
superfiuid density depends on field. Instead it is an in- 
dication that the theoretical model for n(B) does not 
contain all of the relevant physics. On the other hand, 
a "true" measure of the magnetic penetration depth in a 
/iSR experiment has been demonstrated to be given by 
the H — ► extrapolated value of A a b [4l|, [42| • In other 
words, when an inadequate model for n(B) is used to fit 
the TF-/J.SR signal, the difference between A a & and the 
magnetic penetration depth grows with increasing mag- 
netic field. This has been nicely demonstrated in several 
theoretical works [37l. l39l. IH E | . 

The effective length scale X a b measured by [iSR in the 
vortex state of single-crystal superconductors usually ex- 
hibits a strong dependence on magnetic field at low tem- 
peratures. This was first observed in YBa2Cu3 06.95 sin- 
gle crystals [H^], where it was attributed to effects asso- 
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dated with the high anisotropy of the d x i _ y 2-wave su- 
perconducting energy gap A(k). There are two primary 
effects of strong gap and/or Fermi surface anistropy on 

• A nonlinear supercurrent response to the applied 
field H results from a quasi-classical 'Doppler' shift 
of the quasiparticlc energy spectrum by the flow of 
superfluid around a vortex [91| ■ When the Doppler 
shift exceeds the energy gap, Cooper pairs are bro- 
ken, and A increases. Nonlinear effects are par- 
ticularly strong in a d rr 2_ y 2-wave superconductor, 
where even at low T and weak H pair breaking oc- 
curs from the excitation of quasiparticles at the gap 
nodes. A sizeable nonlinear response can also oc- 
cur in s-wave superconductors that have a highly- 
anisotropic energy gap, a smaller energy gap on one 
of the Fermi sheets {i.e. a multi-band superconduc- 
tor), and/or a highly-anisotropic Fermi surface. 

• A nonlocal supercurrent response to the applied 
field H occurs in the limit A << £o> meaning that 
the supercurrent density j(r) depends on the mag- 
netic vector potential A(r) within a volume of ra- 
dius ~ £o around r. Since £o = ^Wp/ttA, the 
value of the superconducting coherence length is 
k-dependent, i.e. it is dependent on both the en- 
ergy gap A(k) and the Fermi velocity i>F(k). In a 
d 2 ,2_ y 2-wave superconductor, nonlocal effects arise 
from the divergence of £o(k) at the gap nodes, but 
they may also occur in isotropic-gapped supercon- 
ductors having strong Fermi surface anisotropy. In 
either situation, the nonlocal response vanishes far 
from the vortex cores, so that only the spatial dis- 
tribution of magnetic field in and around the vortex 
cores is influenced. 

In 1999, a high-field TF-^SR study of YBa 2 Cu 3 6 .95 
[3^ showed that X a b{H) has a sublinear dependence on 
H, in apparent agreement with theoretical predictions 
for nonlocal and nonlinear effects in the vortex state of a 
<i 2 .2_j / 2-wave superconductor 37]. In this case the domi- 
nant contribution to the field dependence of X a b(H) are 
nonlocal effects. With increasing H, the increased over- 
lap of the spatial regions around the vortex cores char- 
acterized by a nonlocal response, reduces the width of 
the fiSK line shape. Since the phenomenological GL and 
London models for B(r) do not account for this, the fit- 
ted value of X a b exceeds the actual magnetic penetration 
depth. 

There have been several /iSR studies on systems that 
exhibit field-induced hexagonal-to-square FLL transfor- 
mations, where the /iSR signals have been fit assuming a 
phenomenological London model developed by Kogan et 
al. 0, loll ]. The Kogan model includes nonlocal correc- 
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tions that stem from Fermi surface anisotropy. As shown 
in figure 1 1 31 the Kogan model adequately accounts for 



FIG. 13: Magnetic field dependence of the fitting parameter 
X a b plotted versus reduced field for LuNi 2 B 2 C [55(, VsSi [641 ]. 
NbSe 2 EH and YBa 2 Cu 3 6 .95 [11]. The reduced field for 
YBa 2Cu3 Q6.95 assumes H C 2 = 1100 kOe, as determined in 
Ref. [l2(|- The results for V 3 Si and LuNi 2 B 2 C come from fits 
to the 'Kogan model' [60l . [6lJ that includes nonlocal effects, 
whereas the data for NbSe 2 and YBa 2 Cu30e.95 come from 
fits assuming an analytical GL model [23( ] . 



nonlocal effects in VsSi at low field, as the fitted X a b 
is field-independent. The same model has also been ap- 
plied to /zSR measurements on borocarbide superconduc- 
tors [HI, [Hf! and NbaSn [g^]. However, in these cases X a b 
still exhibits a field dependence. This suggests that there 
is a highly anisotropic superconducting gap (or multiple 
gaps) that causes the field dependence of X a b — which in 
turn indicates that the theoretical model fails to accu- 
rately portray the field distribution of the FLL. 

Kadono [l2l| has argued that the field dependence of 
X a b measured by /iSR in any type-II superconductor is in- 
dicative of the degree of anisotropy of the superconduct- 
ing order parameter. This seems to be true, although 
Kib(H) is partly influenced by changes in the field decay 
outside the vortex core arising from an increased overlap 
of the quasiparticle core states of neighbouring vortices. 
Some examples of the field dependence of A Q & obtained 
in ^iSR studies are shown in figure [121 For the case of 
VsSi, where the quasiparticles are tightly bound to the 
vortex cores at low reduced field H/H c2 , X a b is at most 
weakly dependent on field. The Fermi surface anisotropy 
responsible for the square FLL that occurs at higher fields 
does not affect X a b(H) at low field. On the other hand, 
the rapid derealization of quasiparticles in NbSe2 with 
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increasing magnetic field appears to influence \ a b(H), as 
it exhibits a stronger field dependence than what is ob- 
served in hi ghly anisotropic YBa 2 Cu 3 06.95. Recently, 
Laiho et al. 44| have shown that the field dependences 



of X a b in VaSi and NbSe2 vanish if an appropriate tem- 
perature and field-dependent cutoff function numerically 
calculated from the quasiclassical Eilenberger equations 
for an s-wave superconductor is used to fit the /iSR mea- 
surements. This is because the low-energy excitations 
of the vortex cores that affect the field distribution are 
handled by the quasiclassical Eilenberger theory. While 
the magnetic penetration depth of an isotropic s-wave su- 
perconductor in the vortex state has no appreciable field 
dependence in the microscopic theory, the vortex core 
size has a strong field dependence due to the intervortex 
transfer of quasiparticles 81 , 8^] • 



MAGNETISM IN HIGH-T C SUPERCONDUCTORS 

At the turn of the 21 st century, neutron scattering, 
NMR and /iSR experiments on high-T c superconductors 
revealed that the application of a magnetic field may 
induce or enhance antiferromagnetic (AF) spin correla- 
tions. In many cases, the AF spin correlations are most 
pronounced in and around the vortex cores where super- 
conductivity is suppressed. These findings suggest there 
is a competing magnetic order that coexists with super- 
conductivity. However, the theoretical implications of 
this hinge on some important details about the experi- 
ments, such as whether the observed magnetism is static 
(or quasistatic) or dynamic, and whether or not static (or 
quasistatic) magnetism is also present in zero magnetic 
field. 



Magnetism in zero applied magnetic field 

High-T c superconductivity emerges from the gradual 
destruction of the insulating AF parent compound by the 
doping of charge carriers into the CuC>2 layers. Shortly 
after their discovery, zero-field (ZF) /iSR experiments 
on high-T c superconductors showed that static electronic 
moments are still present in lightly-doped superconduct- 
ing samples at low temperatures. In particular, an early 
ZF-/zSR study of YBa 2 Cu 3 O y at T = 90 mK showed that 
static magnetism was still present in samples wit h ox - 
gen content y<6.51, but not in a y = 6.54 sample 11221 . 
Similarly, an early ZF-/1SR study of La2- x Sr x Cu04 [1231 ] 
reported static internal fields at low temperature for 
x<0.15. While these experiments suggested that static 
magnetism and superconductivity coexist, there was con- 
siderable worry that the samples studied contained an in- 
homogeneous concentration of cop per m agnetic moments 
due to nonuniform doping of holes 124j . Subsequent ZF- 
/iSR experiments on pure L^-xSr^CuC^ indicated that 
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FIG. 14: Time evolution of the muon spin polarization in 
YBa2Cu3 06.95 at T = 55 K and H = 0, measured w ith t he ini- 
tial polarization P(0) perpendicular to the c-axis |l27l ]. The 
solid curve is the static Gaussian Kubo-Toyabe function with 
A = 0.1144 /is -1 . The inset shows what this function looks 
like beyond the time range of the fiSK signal. Note, the 
data shown here was obtained using a continuous-wave muon 
beam at TRIUMF (Vancouver, Canada). With a pulsed muon 
source, such as that at ISIS (Oxford, United Kingdom), it is 
very feasible to extend this /^SR spectrum up to 20 /is. 



static electronic moments are not pres ent in samples with 
strontium content as high as x = 0.15 [l2 l [l25 l 1l2^. al- 
though the extrapolated T — > K critical value of x for 
the onset of static magnetism has never been accurately 
determined. At dopings above where static magnetism 
is not observed in ZF-/iSR experiments, magnetic fluc- 
tuations persist that are visible by NMR and inelastic 
neutron scattering. 

Whether the frozen spins observed by ZF-//SR in un- 
derdoped superconducting samples coexist with super- 
conductivity or are simply due to phase segregated hole- 
poor regions, has been a central issue of controversy 
in the field. Establishing coexistence experimentally re- 
quires showing that a given sample exhibits both mag- 
netism and superconductivity throughout 100 % of its 
volume. Since magnetic and non-magnetic regions in the 
same sample give distinct /iSR signals with amplitudes 
proportional to the volume of the sample occupied by 
each phase, in principle //SR can determine whether mag- 
netism is present in the entire sample. One constraint, 
however, is imposed by the sensitivity of /iSR to nuclear 
dipole moments. Randomly orientated nuclear dipoles 
cause a Gaussian relaxation of the ZF/iSR time signal, 
as shown in figure [TJ] The solid curve is a fit to a static 
Gaussian Kubo-Toyabe function 



Gf T (t) = | + |(1- A 2 t 2 ) exp^A 2 * 2 ) , (22) 

where A/7^ is the width of the Gaussian field distribu- 
tion at the muon site due to the dense nuclear dipoles 
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(recall 7^ = 0.0852 /is"^ 1 ). Typically, A - 0.1 /is" 1 , 
which corresponds to a field of ^ 1 G at the muon site. 
A muon sitting 20 A away from a copper atom with a 
magnetic moment of ~ 1 /ib will sense a ~ 1 G field, 
whereas a muon sitting further away will not detect a 
field larger than that from the nuclear dipoles. Thus 
if all muons implanted in the sample experience a lo- 
cal magnetic field > 1 G, the nonmagnetic hole-rich re- 
gions cannot be larger than ~ 20 A. T his conclusion was 
reached by Niedermayer et al. [l25l | in ZF-/iSR stud- 
ies of the La2- a; Sr ;c Cu04 and Yi_ a; Ca 3 ;Ba2Cu306 sys- 
tems, where the hole doping in the CuC-2 laye rs is con- 
trolled by cation substitution. Kanigel et al. [l28l ] also 
concluded from ZF-/iSR measurements that phase sep- 
aration in (Ca x Lai_ x )(Bai.75_ ;E Lao.25+x)Cu30j / is on a 
microscopic scale. The maximum size of the nonmag- 
netic hole-rich regions inferred from these experiments is 
nearly equivalent to the minimum possible spatial extent 
of the superconducting regions. The latter is set by the 
superconducting coherence length £0, which is typically 
15-20 A in high-T c superconductors. Thus, the length 
scale of these ZF-//SR experiments is almost, but not 
quite small enough to determine whether magnetism and 
superconductivity coexist in the same sp atial regio n. 
While the ZF-^iSR studies of Refs. [HI, EH pro- 



vided strong evidence for the coexistence of magnetism 
and superconductivity on a nano-scale level, a couple of 
questions about these works have lingered. The first is 
whether superconductivity occurs in the entire bulk of 
the sample. The second question is whether the inher- 
ent cation di sorder th at is present in the systems stud- 
ied in Refs. [l25l Il28l | is responsible for the occurrence 
of electronic moments in superconducting samples. To 
address these questions, more recently, ZF-^iSR, TF- 
fj.SH and bulk magnetization measurements were per- 
formed on both hig h-qua lity polycrystalline samples [1291 ] 
and single crystals [130l | of YBagCu^Ov. In contrast to 
the systems studied in Refs. [125l Il28j . hole doping in 
YBa2Cu30 y is done by varying the oxygen concentra- 
tion in the CuO chain layers, rather than by cation sub- 
stitution. While superconductivity is widely believed to 
take place within the Cu02 planes, the cha i ns th emselves 
are conductive and superconducting 13lL 1 132 1 , due t o 
strong coupling with the nearby CuC-2 planes [l33l . Il34j . 
Furthermore, charge ordering in the CuO chain layers 
indu ces charge-density modulations in the Cu02 layers 
[l35l |. Precision ZF-/zSR measurements on YBa 2 Cu30 y 
show that the c harg e ordering affects the relaxation of 
the [iSR signal [1361 ] . However, this could be due to a 
change in the nuclear dipole field sensed by the muon, 
rather than the formation of electronic moments. More 
importantly, the effects of charge ordering on the ZF- 
/iSR spectrum weakens with decreasing oxygen content, 
making it easy to identify the onset of static electronic 
moments in lightly doped samples. 

The more recent ZF-//SR of studies of high-quality 



YBa2Cu30j, samples [l29l Il30j j are consistent with the 
earlier experiments, in that static electronic moments oc- 
curring in superconducting samples with a hole-doping 
concentration less than p « 0.08 are sensed by all im- 
planted muons. In samples where T c is greater than the 
onset temperature for static electronic moments, TF-/iSR 
measuremen ts ab ove the spin freezing temperature Tf (or 
Tm in Ref. [l3Clj ] ) show that a vortex lattice is formed 
throughout the bulk of the sample. In other words, the 
entire sample exhibits superconductivity. Below Tf, the 
same conclusion cannot be reached, because the mag- 
netism rapidly depolarizes the TF-/xSR signal. However, 
bulk magnetization measurements show no reduction in 
the superconductin g^ res p onse below Tf. Thus, the ex- 
periments of Refs. [l29L Il30l ] appear to establish that 
static magnetism and superconductivity coexist on a mi- 
croscopic scale. Saying any more than this goes beyond 
the spatial resolution of the ZF-/zSR experiments. 

The arrangement of the frozen Cu spins in the coexis- 
tence phase inferred from the ZF-/iSR experiments may 
be compared to that indicated by neutron scattering ex- 
periments. The ZF-/xSR experiments on L^-zSr^CuO^ 
[125j show that a cluster spin-glass (CSG) phase persists 
in the superconducting state and is characterized by an 
onset temperature that decreases continuously across the 
insulator-superconducting boundary. The CSG phase, 
which has also been detected b y 13 ^ La nuclear quadrupole 
resonance (NQR) EM Eli], consists of AF spin 

correlations within domains that are separated by walls 
of hole-rich material, with the easy axis of the stag- 
gered magnetization uncorrelated between clusters. Neu- 
tron scattering experiments on L^-zSr^CuO^ show that 
both static and dynamic incommensurate spin correla- 
tions exist across the insulator-superconducting bound- 
ary |l4dj . The source of the low-temperature static in- 
commensurate spin correlations that give rise to an elas- 
tic neutron scattering peak shifted with respect t o th e 



AF position for 0.2 < x < 0.12 samples [141|, [142J, [14 



is likely the same static (or quasistatic) magnetism de- 
tected by the low-frequency /iSR and NQR measure- 
ments. In fact, Wakimoto et al. |l43| showed that the 
ordered magnetic moment decreases continuously across 
the insulator-supe rcond ucting boundary. In recent years, 
spin-spiral models [l44| have emerged to explain the mag- 
netic incommensurability of La2-a;Sr x Cu04. Liischer et 
al. [l45j have proposed that the CSG phase observed 
at low T by ZF-/iSR results from the freezing of incom- 
patible spin-spiral configurations in neighbouring C11O2 
layers. 

The ZF-/1SR signals for p < 0.08 superconducting 
samples of YBa2Cu30j, at low temperature are charac- 
terized by a r apidly damped oscillation at early times 

[M El Eli 

An example of this is shown in figure [16] 
for YBa 2 Cu306.37- The rapidly damped oscillation is in- 
dicative of very short-range magnetic order. Sanna et al. 
[l29j have argued that this is consistent with nanoscale 
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FIG. 15: Magnetic phase diagrams deduced from ZF-/1SR studies of YBa2Cu30 !/ [1291 ] and La2-iSr x Cu04 [125]. In 
YBa2Cu30j,, static short-range AF correlations are detected below the spin-freezing temperature Tf. In L^-^Sr^CuOzt, a 
'cluster' spin-glass (SG) state appears below T g . In both materials there is a microscopic coexistence of superconductivity (SC) 
and frozen AF correlations in underdoped samples. The insets in each panel show the T— >0 K extrapolation of the transition 
temperatures for static magnetism plotted on a linear scale. 




FIG. 16: Temperature dependence of the ZF-/iSR time signal for YBa2Cu30 H singles crystals recorded wit h the initial muon- 
spin polarization P(0) parallel to the c-axis. The spectra for y = 6. 37 come from the measurements of Ref. [l30t ] . whereas the 
spectra for y = 6.50 and y = 6.57 come from the work of Ref. 158]. The oscillations observed at low T for y = 6.37 indicate 
short-range magnetic order. The ZF-/1SR signals of y = 6.50 and y = 6.57 are identical, temperature-independent, and described 
by relaxation due to nuclear dipole moments and an additional weak exponential relaxation of unknown origin 127]. 



coexistence of stripe-like AF domains and superconduct- 
ing material. The situation is in contrast to the coexis- 
tence of CSG magnetism and superconductivity observed 
in cation-substituted cuprates, which could be argued to 
be an extrinsic property arising from inherent disorder 



Recen t neutron scattering measurements by Stock 
147j on YBa2Cu3C>6.353 single crystals support 



[146 

et al. 

the occurrence of a single phase of coexisting supercon- 
ductivity and short-range spin correlations. Fluctuating 
short-range commensurate spin correlations that develop 
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FIG. 17: fiSK line shapes for I^-^Sr^CuCU single crystals at H = 5 kOe and T = 5 K. For comparison, all line shapes have been 
normalized to their respective peak amplitude n p ^(B). The dashed vertical line indicates the applied field Bq. The background 
signal from muons missing the sample and stopping elsewhere is very small, and hence there is no visible peak at B = Bo. 
The two panels on the right are blowups of the lower portion of the pSR line shapes, with the widths of the line shapes made 
equivalent to account for changes in X a b- The shape of n(B) clearly changes between 0.145 < x < 0.166. 



at temperatures well above T c gradually slow down to 
form a spin glass-like state at low temperature. How- 
ever, in contrast to the ZF-/iSR measurements, no tran- 
sition to magnetic order of any kind is seen in the neu- 
tron measurements. Stock et al. have suggested that 
the magnetic freezing transition t emperatu re observed in 
the ZF-/iSR experiments of Refs. 129l 130l ] is simply the 
temperature at which the slow spin dynamics observed 
by neutron scattering comes into the /iSR time window. 
While this explanation is reasonable, it does not account 
for the oscillation observed in the ZF-/iSR signal at low 
T. A combined /iSR and neutron scattering study on 
the same severely underdoped sample may reconcile this 
discrepancy. 



Field-induced/enhanced magnetism 

NMR ex periments on near-optimally doped 
YBasCusOy Hi5 fl49l HHoj . YBa 2 Cu 4 8 [l5l| and 
Tl 2 Ba2Cu06+5 [l52l ] have detected AF spin fluctu- 
ations in the vicinity of the vortex cores. Likewise, 
low-energy field-induced spin fluctuations have been 
detected by inelastic neutron scattering experiments on 



Lai.837Sr .i63CuO4 |153l |. As the fast spin fluctuation 
rates detected in these experiments fall outside the /iSR 
time window, the field-induced spin dynamics in the 
high-doping regime cannot be studied by /LtSR. On the 
other hand, spin fluctuations slow down considerably 
in the low-doped regime and it is here where further 
information has recently been obtained by //SR. 



La,2- x Sr x CuOi 

There have been a couple of TF-/tSR studies pro- 
viding evidence for frozen spins i n th e vortex cores of 
La2-;cSr x Cu04. 



Kadono et al. |154j investigated the 
temperature dependence of the vortex core size at H = 
2 kOe in both underdoped and overdoped L^-zSr^CuCU 
single crystals. While the temperature dependence of the 
vortex core size at 2 = 0.19 was found to be in qualitative 
agreement with the Kramer-Pesch effect, in the £ = 0.13 
and a; = 0.15 samples the vortex core size increased with 
decreasing temperature and to a value significantly larger 
than expected from measurements of H C 2. The large size 
of the vortex cores in the lower doped samples was argued 
to result from a suppression of superconductivity around 
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FIG. 18: Temperature dependence of the exponen- 
tial relaxation rate A from fits of the TF-/1SR sig- 
nal for Lai.855Srn.i45Cu04 at H = 5 kOe (circles) and 
Lai.8 24 Srp. 176 CuQ4 at H = 5 kOe (triangles) to equation (123|) 
[l58| . Temperature dependence of the exponential relaxation 
rate Ai (squares) from fits of the 'envelope' of the TF-/^SR 
signa l for Lai.8sSro.i2Cu 04 a t H = 2 kOe, 10 kOe and 50 kOe 
[l6ll ]. Note, A from Ref. [I58j | is associated with magnetism at 
the centre of the vortex core (r = 0), where superconductivity 
is fully suppressed. 



the vortex due to the formation of AF spin correlations. 
At the time, this seemed co mpa t ible with neutron exper- 



iments on La2- :c Sr ;c Cu04 [1551 . 1 1 561 ] showing enhanced 



static AF spin correlations in underdoped samples. How- 
ever, Kadono et al. speculated that the spin correlations 
responsible for the large vortex core size at low T and low 
H in La2- 2: Sr :E Cu04 with x = 0.13 and a; = 0.15 are not 
static, but dynamic. This interpretation was supported 
by a subsequent neut ron study of an x = 0.144 sample 
by Khaykovich et al. 1571 ] that revealed static magnetic 



order only above H w 30 kOe. However, not ruled out 
by these studies is the possibility that there is disordered 
static magnetism at low fields. 

The latter possibility was considered in a more re- 
cent /iSR study of La 2 - x Sr :E Cu04 single cryst als that 



are free of static magnetism in zero field |158j . Visu- 
ally, evidence for static magnetism confined to the vor- 
tex cores can be seen in the TF-/xSR line shape as an 
unusual high-field 'tail' and the possible appearance of 
a low-field tail (depending on the magnitude and ori- 
entation of the fields in the vortex core), but with no 
other change in the shape of n(B). As shown in fig- 
ure [TTl the high-field tail of the /iSR line shape for 
La 2 - 2; Sr :E Cu04 changes at a Sr doping 0.145 <x< 0.166. 
However, similar changes of the /iSR line shape could 
also arise from a 3D-to-2D crossover of the FLL. For ex- 
ample, in highly anisotropic Bi2+a;Sr2- x CaCu208+a Lee 
et al. 159] showed that above a crossover field, random 



pinning-induced misalignment of 2D 'pancake' vortices 
in adjacent layers narrows and reduces the asymmetry of 
the /xSR line shape. For La2_2.Sra.CuO4, a similar sce- 
nario would imply that a reduction in Sr substitution 
significantly softens the flux lines and induces a degree 
of disorder sufficient to misalign the 2D vortices in ad- 
jacent layers. However, a combined /xSR and neutron 
study of Lai.gnSro.irjCuCU by Divakar et al. [z2] demon- 
strated that the vortices are 3D-like even at this lower Sr 
concentration. This conclusion is supported by neutron 
scat terin g experiments on Lai.gnSrn.ioCuC^ by Lake et 
al. [l60l | . which show that the field-induced static mag- 
netic order is 3D. 

While the onset of magnetism is a more likely cause of 
the change in the /xSR line shape of La2-. r Sr r C uC>4 with 
doping, the neutron study of Khaykovich et al. [1571 ] im- 
plies that the line sha pe fo r x = 0.145 measured at low 
T and low H in Ref. 1581 ] does not result from static 



magnetic order. However, the TF-/1SR precession signals 
were found to be well described by the following polariza- 
tion function, which allows for spin-glass like magnetism 
in and around the vortex cores 



P(t) 



exp( 



-Ae 



-{n/uy 



i)cos[7„fl(ri)«]. (23) 



Here the sum is over all sites in the real-space unit cell of 
the FLL and 5(r,) is the local field at position n = (xi,yi) 
with respect to the vortex center. The distribution of 
fields at each site is assumed to be Lorentzian, and the 
relaxation rate A due to the magnetism is assumed to fall 
off as a function of radial distance from the vortex core 
centre on the scale of £ a {,. The temperature dependences 
of A for La2_ a; Sr :I ;Cu04 with x — 0.145 and x = 0.16 at 
H = 5 kOe are shown in figure [TH] The increase of A 
with decreasing temperature observed for the x = 0.145 
sample indicates a slowing down of spin fluctuations in 
the vortex-core region. 

It is important to emphasize that these fits to the TF- 
/iSR time spectra don't in themselves prove the existence 
of spin-glass order in the vortex cores. Other models 
based on different microscopic physics, including differ- 
ent forms of magnetism, may give equally good or better 
fits to the time-domain signals. The conclusion that mag- 
netism occurs in the vortex cores really comes from the 
changes in the TF-/1SR line shape that are visually appar- 
ent as a function of temperature, magnetic field and hole 
doping. When the changes as a function of these three 
independent variables are considered together, there are 
really only two logical interpretations of the data. As al- 
ready mentioned, experiments that have probed the di- 
mensionality of the vortices in La2_ a; Sr :I ;Cu04 more or 
less rule out the cause being a vortex transition. In 
the next section, experimental evidence that also argues 
against this being the interpretation for similar changes 
of the TF-/xSR line shape of YBa2Cu30 a is discussed. 

With this in mind, the results of the fits for 
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FIG. 19: jUSR line shapes for YBa2Cu30 !/ single crystals at H = 5 kOe and T — 2.5 K. For comparison, all line shapes have 
been normalized to their respective peak amplitude n p ^(B). The dashed vertical line indicates the applied field Bo, and the 
bump or peak near this line (most visible for y = 6.50) is the background signal. The two panels on the right show blowups 
of the lower portion of the /iSR line shapes, with the widths of the line shapes made equivalent to account for changes in X a b- 
The shape of n(B) clearly changes between 6.50 < y < 6.57. 



Lai.855Sro.i45Cu04 may be compared to recent TF -^SR 



measurements of Lai.8 8 Sro.i2Cu04 by Savici et al. [1611 ] . 
showing a field-induced enhancement of static or qua- 
sistatic magnetism. In contrast to the TF-^SR stud- 
ies discussed so far, information on field-induced mag- 
netism was obtained in this work by fitting the 'enve- 
lope' of the muon spin precession signal at early times to 
a two-component exponential function E\ exp(— K\t) + 
E% exp(— A 2 t), where the first and second terms describe 
contributions to the [iSR signal from volume fractions 
of the sample with and without quasistatic magnetism, 
respectively. This was possible for Lai.8sSro.i2Cu04, be- 
cause the field-enhanced magnetism occupies a volume 
of the sample much larger than the volume of the vor- 
tex cores. In particular, Ei/E 2 > 60 % co mpared to 
Ei/E 2 <l % Lai. 855 Sro.i4 5 Cu04 in Ref. [l58| . As shown 
in figure [THJ the temperature dependence of the relax- 
ation rate A of Lai.855Sro.i45Gu04 determined from fits 
of the full TF-/xSR precession signal with equation (|23j) is 
similar to the temperature dependence of the relaxation 
rate Ai determined from fits to the envelope function of 
Lai.88Sro.i2CuC>4 at low H. This seems to suggest that 
the static (or quasistatic) magnetism that occurs in a 
large volume of Lai.8 8 Sro.i2Cu04 in a magnetic field is 



similar to the magnetism found near the vortex cores in 
Lai.855Sro.i45Cu04. However, this is unlikely the case 
at low temperatures. In contrast to Lai.sssSro.usCuC^, 
static magnetic c orre lations occur in Lai.8 8 Sro.i2Cu04 
even in zero field |l56| . Thus the fast relaxation rate ob- 
served in Lai.ssSro.^CuCU in the vortex state at low T 
can be attributed to static magnetic order. 



The study by Savici et al. also showed that quasistatic 
magnetism persists to temperatures well above T c at high 
H . This implies that vortices are not required to nucleate 
static magnetism in superconducting samples. Instead 
the vortices act to suppress spin fluctuations and en- 
hance the magnetic volu me fraction. The TF-/iSR exper- 
iments by Savici et al. 161] provide strong evidence for 
these effects, as well as neutron scatte ring experiments on 
Lai.8 8 Sro.i2Cu04 by Katano et al. 156] showing field- 
enhanced AF spin correlations at H — 100 kOe and low T, 
and Raman scatteri ng ex periments on Lai.8 8 Sro.i2Cu04 



by Machtoub et al. 1621 ] that can be explained in terms 



of a field-enhanced AF volume fraction. 
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YBa 2 Cu 3 O y . 

Miller et al. |l63| reported evidence from TF-^tSR 
measurements on YBa 2 Cu306.5o at H > 10 kOe for 
the occurrence of sta tic A F order in the vortex cores. 
While Kadono et al. (Ull were able to fit low field TF- 
/iSR spectra for underdoped La2- K Sr :E Cu04 assuming 
the spatial field profile B(r) of the FLL is adequately 
described by the London model with a Lorentzian cutoff 
factor F(G) = exp(~^/2G£_), fits of the ^SR time spec- 
tra for YBa2Cu30g.5o assuming a conventional Lo ndon 
or GL model were found not to converge. In Ref. [l63j 
the TF-^SR signal from YBa2Cu306.so could be fit as- 
suming static AF order in the vortex cores commensurate 
with the underlying crystal lattice, even though to date 
field-induced static magnetic order has not been detected 
in YBa2Cu30j, by neutron scattering. The magnitude of 
the staggered field sensed at the muon stopping site(s) 
was determined to be ~ 18 G at low T. However, this 
analysis still required the value of £ a b to be constrained, 
which suggests an insensitivity to the high-field tail of the 
/iSR line shape or that the assumed theoretical model is 
inappropriate. 

Recently, the doping, temperature and magnetic field 
dependences of the /iSR line shape for YBa^ CusO^ sin- 
gle crystals were studied in some detail jl58| . fiSR line 
shapes for YBa2Cu30g.46 and YBa2Cu30g.5o were found 
to be compatible with spin-glass magnetism appearing in 
and around the vortex cores (some of these measurements 
are shown in figure I19p . This interpretation of the /iSR 
signal explains why frozen spins have not been detected 
in the vortex cores of YB^CvL^Oy by neutron scatter- 
ing, which is insensitive to static magnetism with random 
spatial correlations. As in the case of La2- x Sr x Cu04, 
an explanation of the observed changes in the TF-/xSR 
line shape in terms of a crossover in the dimensionality 
of the vortices is unlikely. Josephson plasma resonance 
measurements by Dulic et al. [801 ] indicate that the vor- 
tices in YBa 2 Cu306.5o are 3D-like at low temperatures. 
Consistent with this finding, mutual inductance measure- 
ments show that even severely underdoped YBa2Cu30 y 
becomes quasi-2D only near T c [30l ]. 



Pr2- x Ce x Cu04 and Pr\- x LaCe x Cu04. 

The first reported experimental evidence for field- 
induced magnetic order in an electron-doped cuprate su- 
perconductor came from a TF-/xSR study of underdoped 
P^-zCezCuO^ single crystals 164J. In zero field, ran- 
dom magnetism was detected, but a remarkably weak 
applied field of only 90 Oe resulted in the occurrence 
of static AF order thoughout the volume of the sample. 
The onset of bulk AF order below T c is visually apparent 
from the observed positive field-shift of the entire /iSR 
line shape (see figure |2"0")) . In this case, the entire line 
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FIG. 20: Temperature dependence of the fiSK line shape in 
an u nderdoped P^-xCe^CuC^ single crystal at H = 90 Oe 
[l64l |. Inset: Temperature dependence of the average internal 
magnetic field Bo. The blue and red colour scheme denotes 
measurements taken above and below T c , respectively. 



shape shifts to higher B because the AF order enhances 
the local magnetic field at all of the 0(3) muon stopping 
sites. 

The occurrence of weak field-induced AF or- 
der was subsequently identified in superconducting 
Pr 1 _ x LaCe 2: Cu04 with x = 0.11 and x = 0.15, in TF-/xSR 
studies by Kadono et al. [166, 167]. In these experi- 
ments it was argued that the muon senses the antifer- 
romagnetism indirectly through a transferred hyperfine 
coupling between the fi + spin and the Cu moments via 
the Pr ions, rather than by direct detection of the dipolar 
fields o f the Cu m oments at the muon site as suggested 
in Ref. [l6ill65j l. 

Single crystals of Pri_ x LaCe x Cu04 can be made much 
larger than single crystals of Pr2_ :c Ce :r Cu04, and this 
has allowed neutron scattering studies of this com- 
pound. Fujita et al. used neutrons to study the same 
Pri-^LaCe^CuO ^ sa mples that Kadono et al. investi- 
gated with //SR |l68| . In zero field, weak AF order was 
detected in the x = 0.11 sample, in contrast to ZF-/iSR 
measurements that fo und only random magnetism due 
to small Pr moments [l69l |. This apparent discrepancy 
suggests that the AF correlations in zero field are fluctu- 
ating on the /xSR time scale. The neutron experiments 
by Fujita et al. also demonstrated a field-induced en- 
hancement of AF order in the x = 0.11 sample, whereas 
static AF order was not observed at all in the x = 0.15 
sample. The latter result contrasts the TF-/zSR experi- 
ments of Kadono et al. who reported field-induced AF 
order in the x = 0.15 sample, but with a reduced ordered 
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FIG. 21: Schematic of theoretically predicted zero-temperat ure p hase diagrams for (a) 2D vortices (adapted from figure 1 of 
Ref. |l7l| ). and (b) 3D vortices (adapted from figure 1 of Ref. [173| V Here a is a control parameter (e.g. carrier concentration). 

and <j> denote the expectation values of the superconducting and competing order parameters, respectively. The solid blue 
curve represents a quantum phase transition (QPT) and the solid blue dot denoted ct c is a quantum critical point (QCP). Note 
that the QPT for the 2D vortex case is a crossover within the coexistence phase of the 3D vortex phase diagram (red dashed 
curve). Below the crossover curve, the competing order is spatially inhomogeneous. In addition, the zero field extrapolation of 
the QPT in (b) is an 'avoided' QCP (open blue circle). 



local moment. More recently, Kang et al. have detected 
static SDW order and residual AF order in underdoped 
Pri_ a; LaCe a ;Cu04 at zero field by neutron sc atteri ng, but 
not in optimally doped Pri_2.LaCej.CuO4 [170|. Fur- 
thermore, they showed that a c-axis aligned magnetic 
field enhances the SDW order, but not the AF order 
in underdoped Pr 1 _ a; LaCe ; ;;Cu04, and does not induce 
static magnetic order of either kind in optimally doped 
Pri_ x LaCe a: Cu04. 



Piecing together the information from /iSR 

ZF-/iSR and TF-,uSR studies of cuprates have es- 
tablished that superconductivity and static (quasistatic) 
magnetism coexist on the low-doping side of the super- 
conducting 'dome'. In quality samples, this coexistence 
is on a nanometer scale at zero magnetic field. Field- 
induced static magnetism has been detected at low tem- 
peratures in samples where the spin fluctuation rate at 
H = falls outside the /zSR time window. In these 
samples, disordered static magnetism is induced by weak 
magnetic fields. 

As discussed in Ref. [l4dj j . neutron scattering studies 
of La2_ K Sr a; Cu04 and La2Cu04+j, support a theoreti- 
cal model of competing superconducting and mag netic 



order parameters proposed by Demler et al. [17l|. In 
this model the pure superconductor undergoes a field- 
dependent quantum phase transition (QPT) to a state of 
coexisting superconductivity and magnetic order (see fig- 
ure [5T^,). The competing magnetic order is stabilized at 
the QPT by the suppression of superconductivity associ- 



ated with the formation of vortices. In agreement with 
the theory, the neutron e xper iments have detected field- 
induced magnetic order [1571 ] and have shown that the 
ordered magnetic moment grow s with inc reasing mag- 
netic field exactly as predicted [l55l 174 1. The theo- 
retical model of Demler et al. is not exclusive to the 
cuprates, as it also appears to describe the occurrence of 
a field-induced coexistence phase of magnetic order and 
superconductivity in CcRhIn 5 |l72l |. 

In actuality the theoretical model of Demler et al. is 
somewhat oversimplified. In particular, the model as- 
sumes that the vortices are 2D. As a finite size system, 
an isolated 2D vortex cannot support static magnetic 
order. It is only when the 2D vortices strongly couple 
to their neighbours is static magnetic order stabilized. 
Thus the coexistence phase in the theory of Demler et 
al. is characterized by a competing order parameter that 
is nearly spatially uniform throughout the sample. How- 
ever, as mentioned earler, neut ron s cattering studies of 
Lai.9oSro.ioCu04 by Lake et al. [l6C| show that the mag- 
netic o rder in the vortex state is in fact 3D. As explained 
in Ref. [l60l ] , this implies that there is significant coupling 
between the Cu02 planes, so that the vortices themselves 
are 3D. As discussed earlier, this same conclusion was 
reached in a TF-/iSR study of Lai.goSro.ioCu04 by Di- 
vakar et al. 



72] 



Kivelson et al. [1731 ] have extended the theory of Dem- 
ler et al. to account for the coupling between Cu02 
layers, and showed that in contrast to an isolated 2D 
vortex, a competing order parameter could be stabilized 
in an isolated 3D vortex line. This leads to a modifi- 
cation of the phase diagram proposed by Demler et al., 
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FIG. 22: Generic phase diagram inferred from /iSR studies on high-T c superconductors. The acronyms for the different phases 
are defined as follows: AF = antiferromagnetic phase, SM = static magnetic phase (where magnetism may be disordered), SC 
= superconducting phase, and SC + SM = coexisting superconducting and static magnetic phases. T c (p), Tn(p), Tf(p), H C 2(p) 
and pmic denote the superconducting critical temperature, the Neel temperature, the spin freezing temperature, the upper 
critical magnetic field and the T — >0 metal-to-insulator crossover (MIC), respectively. TF-/iSR measurements indicate that the 
crossover at pmic extends below H C 2(p), separating a superconducting phase with fast spin fluctuations from a superconducting 
phase that coexists with static magnetism. ZF-/1SR measurements indicate that superconductivity and static magnetism coexist 
on a nanometer scale at H — 0, but only up to a carrier concentration p that is less than pmic- Above the dashed red curve, 
the static magnetism is uniformly distributed throughout the sample. 



such that the true QPT is to a phase in which a spa- 
tially inhomogeneous competing order coexists with su- 
perconductivity (see figure [2lb) . The TF-^S R experi- 
ments on La 2 _ I Sr K Cu0 4 (0.12 < x < 0.166) [HI HH 
and YBa 2 Cu 3 O y (y < 6.57) [HI [HI] strongly support 
the existence of this phase. The observed changes in the 
/iSR line shape are consistent with the onset of disor- 
dered static magnetism in and around weakly interact- 
ing vortices. These experiments imply that the com- 
peting order parameter stabilized at the QPT is the 
mean squared local magnetization. With increasing mag- 
netic field, stronger overlap of the magnetism around 



neighbouring vortices may lead to a co-operative bulk 
crossover to long-range magnetic ord er, as is apparently 
the case in Lai.g56Sro.i44Cu04 [l57l |. However, so far 
this has not been observed in YBa2Cu3 0j, where the 
QPT identified by /iSR is closer to the low-doping side of 
the superconducting 'dome'. While the QPT in electron- 
doped cuprates has not been accurately determined by 
/iSR or neutron scattering, experiments to date suggest 
the coexistence phase is dominated by static magnetic or- 
der rather than disordered static magnetism — likely due 
to close proximity of the AF and superconducting phases. 

The TF-/^SR experiments on L^-xSr^CuC^ and 
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YBa 2 Cu 3 0y of Ref. [158J indicate that the QPT in fig- 
ure [2Tb occurs at the critical doping below which the 
temperature dependence of the normal-state resistivity p 
at H > H C 2 and low-T changes from metallic behaviour 
(dp / dT > 0) t o an unusual insula ting behaviour (dp/dT < 

This low-temperature 
crossover has also been de- 



0) [17.. 



i7a [ml . E7I IM Eioj 

metal-to-insulator (MIC) 
tected indirectly by thermal co nductivity me asurements 
as a function of magnetic field |l8ll . Il82t Il83l ] . Although 
the MIC is a generic property of high-T c cuprates, it oc- 
curs at a non-universal carrier concentration. The onset 
of static magnetism below the MIC detected by TF-/zSR 
is compatible with the localization of charge. In zero- 
field, deviations from a linear dependence of the resistiv- 
ity on T are observed in underdoped samples, including 
upturns of the resistivity with decreasing temperature 
in very underdoped sup erconduct ing samples. Theoret- 
ically, Marchetti et al. 184 Il85| have argued that the 



cause of the MIC in both zero field and strong magnetic 
fields is t he sa me, and is due to the onset of short-range 
AF order |184j | . This seems compatible with the ZF-/^SR 
studies that have detected static magnetism in under- 
doped samples. 



A key prediction of the theory of Ref. [173l | is that there 
is an 'avoided' QCP at H — 0. In other words, the QCP 
lies at a lower doping than one expects from the extrapo- 
lated H^O location of the QPT at finite magnetic field. 
In La2- a; Sr a ;Cu04, the spin freezing trans ition tempera- 
ture T„_at H = deduced from ZF-/xSR [H| Elf and 
NQR |l38| experiments extrapolates to zero below the 
critical doping for the MIC at x « 0.16 [HI [lSlL Ell ( see 
figurefTS]). The same is true for YP^CuaOy, where static 
magnetis m associated with Cu spins is not observed by 
ZF-^SR [HI EH, EH, E&1 in samples immediately be- 
low the critical doping for the MIC at j/«6.55 [l83(. 



Figure [22] shows the generic magnetic phase dia- 
gram suggested by ZF-/iSR and TF-/iSR experiments on 
cuprates. In this figure the QCP and 'avoided' QCP are 
shown as solid and open circles, respectively. Note that 
the zero-field QCP here is distinct from the universal 
QCP near p = 0.19 that has been suggested by Tallon 
and others 126l 186l 188l |. In other words, there may in 
fact be two unrelated quantum critical points under the 
superconducting 'dome' that strongly influence the phys- 
ical properties of a large portion of the cuprate phase dia- 
gram. While neutron scattering and /iSR measurements 
both support a generic phase diagram ruled by the close 
competition between superconducting and magnetic or- 
ders, it is still to be determined whether some magnetism 
is also responsible for the 'glue' that binds the supercon- 
ducting carriers. 
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